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Abstract 



Plasma pressure distribution in the inner magneto- 
sphere is one of the key parameters for understand- 
ing the main magnetospheric processes including ge- 
omagnetic storms and substorms. However, the pres- 
sure profiles obtained from in-situ particle measure- 
ments by the high-altitude satellites inside the plasma 
sheet and other regions of the magnetosphere do 
not allow tracking the pressure variations related to 
the storms and substorms, because a time interval 
needed to do this generally exceeds the characteris- 
tic times of them. On contrary, fast movement of 
low-altitude satellites makes it possible to retrieve 
quasi-instantaneous radial or azimuthal profiles of 
plasma pressure along the satellite trajectory, using 
the precipitating particle flux data in the regions of 
isotropic plasma pressure. For this study, we used the 
low-altitude polar-orbiting Aureol-3 satellite data for 
plasma pressure estimation, and the IGRF, Tsyga- 



nenko 2001 and Tsyganenko 2004 storm time geo- 
magnetic field models for the pressure mapping into 
the equatorial plane, and for the evaluation the cor- 
responding volume of the magnetic flux tube. It was 
found that during quiet geomagnetic condition the 
radial pressure profiles obtained coincide with the 
profiles, obtained previously from the high-altitude 
measurements. On the contrary, the plasma pres- 
sure profiles change significantly during the develop- 
ment of storms and substorms. Nevertheless, three 
geomagnetic field models gave significantly different 
geomagnetic field profiles, those points out the neces- 
sity to develop a magnetically self-consistent model 
for description of the inner magnetosphere geomag- 
netic field. However, the common features observed 
for all models are: during geomagnetic storm the 
plasma pressure profiles became sharper; the posi- 
tion of the maximum of plasma pressure corresponds 
to expected one for given Dst minimum; the max- 
imum value of inner magnetosphere static pressure 
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correlates with the solar wind dynamic pressure. In- 
crease in the plasma pressure profiles indicates the 
possibility to consider the interchange instability as 
one of important factors for the development of the 
main phase of geomagnetic storm. 

1 Introduction 

As it is well known, strong geomagnetic storms take 
place when the interplanetary magnetic field (IMF) 
has a prolonged southward orientation. During 
storms a powerful ring current is developed, the in- 
tensity of which is commonly evaluated by Dst index, 
introduced by M. Sugiura (1964). Contribution of the 
ring current into the Dst variation is determined by 
the total inner magnetospheric energy content of the 
trapped particle population, which depends on parti- 
cle supply, energization and losses. The inner magne- 
tospheric energy content is generally evaluated using 
the Dessler-Parker-Sckopke-Tveskoy relation (Dessler 
and Parker, 1959; Sckopke 1966; Parker, 1996; Tver- 
skoy, 1997). Recently, this relationship was modified 
by Antonova (2002) and Liemohn (2003), taking into 
consideration that the value of plasma pressure at 
the external boundary of the ring current region is 
not equal to zero. However, during last decade it 
was proposed that the tail current makes significant 
contribution into Dst variation. This statement was 
based on the fact that the measured by Dst storms 
strength or its rate of strengthening can be decreased 
during substorm expansion phase (Iyemori and Rao, 
1996). It was also found that, on average, during the 
geomagnetic storm main phase, Dst decreases more 
slowly after a substorm expansion phase starts than 
before, and that on average, during the storm's re- 
covery phase, it increases more rapidly (Iyemori and 
Rao, 1996, McPherron, 1997). Nevertheless, the tra- 
ditional interpretation of Dst as a measure of ring cur- 
rent intensity can be reestablished by modifying the 
topology of high latitude current systems (Antonova, 
2004). As it is well known, energetic particles move 
around the Earth near equator in the midnight sec- 
tor and at high latitudes near noon in the region 
of quasi-trapping. Antonova and Ganushkina (2000) 
suggested that the existence of daytime high-latitude 



field minima leads to the splitting of the high latitude 
ring current into two branches. The suggested cur- 
rent system was named the cut-ring current (CRC) 
system. It was shown that the transverse currents, 
created by the gradients of the plasma pressure at 
L more than 7 till the magnetopause, are concen- 
trated far from the equatorial plane in the daytime 
sector of the magnetosphere, being a high latitude 
continuation of ordinary ring current, distorted due 
to daytime compression of the magnetosphere. This 
makes it possible to turn back to the traditional in- 
terpretation of Dst as a result of ring current devel- 
opment, taking into account energy content of the 
region of quasi-trapping in the ring current energy 
budget that includes the nighttime part of the plasma 
sheet till aproximately 12 Re (Radius of the Earth) 
during quiet conditions and its daytime continuation. 
The results of Lyons et al. (2003) made a signifi- 
cant contribution in understanding of processes lead- 
ing to the decrease in the absolute value of Dst or 
in the rate of strengthening of geomagnetic storms. 
In particular, they have shown that substorm injec- 
tions lead to the decrease in the plasma pressure at 
the geocentric distances 10-13 Re. Pressure reduc- 
tions after the substorm expansion phase onset at 
the geostationary orbit were also observed earlier by 
Roux (1985) and Kozelova et al. (1986). So, we can 
assume that plasma transport by the large-scale elec- 
tric field during substorm growth phase leads to the 
increase of Dst. Otherwise, losses of particles after 
the substorm expansion phase onset (probably in the 
tailward direction) explain the decrease of ?Dst? or 
the decrease of the rate of strengthening of the ge- 
omagnetic storm. So, the correct consideration of 
pressure balance in the inner magnetosphere is very 
important for understanding of the ring current dy- 
namics during geomagnetic storms. However, for the 
correct estimation of this balance it is necessary to 
include the ionospheric sources of ions. In partic- 
ular, it was found that the ions of oxigen provide 
more than aproximately 40 percent of particle en- 
ergy density during the main phase of great storms 
and aproximately 20 percent during small to mod- 
erate storms (Daglis, 1997, Pulkkinen et al., 2001). 
Distribution of plasma pressure in the magnetosphere 
of the Earth has been studied extensively during last 
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decades. In particular, Lui et al. (1987), Lui and 
Hamilton (1992) obtained profiles of plasma pressure 
using in situ particle measurements onboard the high- 
altitude AMPTE/CCE satellite. Later, De Miche- 
lis et al. (1999) also obtained bi-dimensional distri- 
bution of plasma pressure in the equatorial plane, 
using the data of the same satellite. Lui (2003) 
studied the magnetic local time (azimuthal) asym- 
metry in the inner magnetosphere plasma pressure. 
Wing and Newell (1998, 2000) reproduced statisti- 
cal bi-dimensional distribution of plasma pressure in 
the equatorial plane at aproxiamtely 10 Rc using 
the low-altitude DMSP scries satellites and deter- 
mined the corresponding Region 1 field-aligned cur- 
rent distribution. Stepanova et al., (2002) modified 
the technique proposed by Wing and Newell (1998, 
2000), considering the presence of a field-aligned po- 
tential drop in the auroral geomagnetic field lines, 
and obtained footprints of plasma sheet pressure pro- 
files in the equatorial plane for very short time in- 
tervals (minutes) for quiet geomagnetic conditions 
and during different substorm phases. This is im- 
possible to do using only the high-altitude satellite 
measurements. Stepanova et al., (2004a) evaluated 
the quasi-instantaneous value of azimuthal plasma 
pressure gradient and corresponding values of field- 
aligned currents, using oblique passages of the low- 
altitude Aureol-3 satellite, that agreed with typical 
values of Region 1 and 2 field-aligned currents. Study 
of plasma pressure profile during magnetic storms is 
constrained because of appearance of energetic parti- 
cle fluxes, mainly relativistic electrons that can cause 
significant distortions in the instrument functionality 
and even the satellite devices. The main contribu- 
tion to the storm time ring current plasma pressure 
is introduced by ions with energies 100 keV, approx- 
imately, (see Lyons and Williams, 1982; Williams, 
1983). However existence of CRC requires analysis 
of particle fluxes with energies of 10 keV, approxi- 
mately, forming part of particle population responsi- 
ble for part of Dst variation. In this paper we use the 
methodology developed by Stepanova et al., (2002) 
for studying the evolution of the radial pressure pro- 
files, obtained by the Aureol-3 satellite before, during 
and after the March 1-8, 1982 geomagnetic storm and 
establish the upper limit of the increase in steepness 



of low energy the part of magnctospheric plasma pres- 
sure profile connected to the ion fluxes with energies 
1-22 keV. We try to show that this limit is determined 
by the development of interchange (flute) instability. 

2 Evolution of the radial 
plasma pressure profiles 
during March 1-8, 1982 
geomagnetic storm 

Despite significant efforts concentrated on the deter- 
mination of distribution of the plasma pressure in 
the inner magnetosphere, the evolution of this dis- 
tribution during geomagnetic storms is not estab- 
lished. As was mentioned before, we propose to ob- 
tain nearly-instantaneous radial plasma pressure pro- 
files at the external part of the ring current using 
low-altitude satellite data. The main difficulties of 
this method consist in the necessity to use some ge- 
omagnetic field model for mapping of the measured 
plasma pressure profile into the equatorial plane and 
possible anisotropy of particle fluxes. High apogee 
satellites obtain the plasma pressure in situ, never- 
theless they move comparatively slowly, and it takes 
hours to obtain a radial pressure profile. This time is 
too large to guarantee what the profile obtained re- 
mains unchanged, and this approach can not be used 
for geomagnetic storm studies. Use data of low or- 
biting satellites makes it possible to restore plasma 
pressure profiles only in the regions of isotropic pres- 
sure. Therefore it can give only partial information 
about the plasma pressure profile. Nevertheless, this 
kind of information gives the possibility to evaluate 
the lower limit of the value of plasma pressure in the 
equatorial plane and steepness of radial plasma pres- 
sure profile. 
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Solar wind parameters for March 1-8, 1982 geomagnetic storm 
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Figure 1: Variations of the solar wind and Dst index between February 25 and March 12, 1982. From top to 
bottom: IMF B x , B y , B z , in GSE, the solar wind velocity and number density, and Dst index. Data were 
downloaded from the OMNI 2 data base ( htt p: / /omniweb.gsfc.nasa.gov/[ ) 



4 




MLT 20.07 20.02 19.97 19.89 19.81 19.69 19.50 

ILAT 51.43 55.22 59.16 63.21 66.87 71.11 75.41 

ALT 557 528 502 478 460 442 427 

L 2.57 3.07 3.80 4.92 6.48 9.54 15.76 

Figure 2: The electron precipitating fluxes (a), the field- aligned current densities and field- aligned potential 
drops (b), the ion precipitating fluxes (c), the ion concentration at the ionospheric altitudes and the ion 
temperature (d), the ion concentration in the magnetosphere and the plasma pressure (e), obtained during 
the March 2, 1982 auroral oval crossing. 
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Here we use the same methodology as in 
(Stepanova et al., 2002) to study the low energy 
plasma pressure profile during the March 1-8, 1982 in- 
tense geomagnetic storm (minimum Dst 211 nT). The 
interplanetary magnetic field components in GSE and 
the solar wind velocity and number density and Dst 
index variations arc shown in Figure 1. Pressure 
profiles were obtained taking into consideration the 
influence of field-aligned potential drop on the ion 
precipitations. This was carried out simultaneously 
with the analysis of the meridional distribution of 
the electron precipitating flux using the polar orbit- 
ing (perigee: 410 km, apogee: 2000 km, inclination: 
82.5), three-axis-stabilized Aureol-3 spacecraft. Par- 
ticle data with a time resolution of 1.6-3.6 s were 
provided by the SPECTRO instruments within the 
0.02-22 keV energy range (Bosqued et al., 1986). 

We analyzed 13 radial profiles of plasma pressure 
obtained during auroral oval crossing in the sector be- 
tween 19 and 23 hours of magnetic local time (MLT) 
before, during and after the March 1-8, 1982 geomag- 
netic storm. The procedure of low energy plasma 
pressure profile retrieval is illustrated in Figure 2. 
In this Figure we can see the electron precipitat- 
ing fluxes (a), the field-aligned current densities and 
field-aligned potential drops (b), the ion precipitat- 
ing fluxes (c), the ion concentration at the ionospheric 
altitudes and the ion temperature (d), the ion con- 
centration in the magnetosphere and the low energy 
plasma pressure (e), obtained during the March 2, 
1982 auroral oval crossing at 20 MLT, approximately. 

This procedure is described in detail in (Stepanova 
et al., 2002, 2004ab). It is based on the as- 
sumption that the ion anisotropy in the equato- 
rial plane is low in the range between 22 and 02 
MLT (De Michelis et al., 1999), and the ion dis- 
tribution functions are nearly Maxwellian. As re- 
ported by Antonova and Tverskoy (1975), the field- 
aligned potential drop decreases the ion concentra- 
tion at the ionospheric altitude and docs not af- 
fect the ion temperature. In this case the ion 
particle flux measured by a low-altitude satellite 
is 1(E) = [n JOIl /(2 1 /2 TO ^3/2Tf/^] Eexvi-E/Ti). 

Here is the ion energy, E is the ion energy, mi is 




the ion mass, and rii, 



, exp (— eV/Ti), where 



{n mag ) and temperature (T^), obtained from precip- 
itating ion fluxes are shown in Fig. 2(d). Figure 
2 (e) shows that low energy plasma pressure has a 
clear maximum 1.2 nPa at 07h32m30s UT. After 
that the pressure decreases till 0.1-0.2 nPa, that sat- 
isfies to the total pressure balance conditions in the 
tail lobes, assuming that in this region the contribu- 
tion of energetic particles in the plasma pressure is 
not significant, and that the value of the magnetic 
field in the lobes is 20 nT, except a small variation 
at 07h34m00s UT. 

The main difficulty in the reconstruction of the 
low energy plasma pressure profile is how to map the 
pressure obtained into the equatorial plane. Figure 3 
shows the geomagnetic field lines, obtained using the 
International Geomagnetic Reference Field (IGRF), 
Tsyganenko 2001 and 2004 geomagnetic field models 
(Macmillan ct al. (2003), Tsyganenko 2002ab, 2005). 
IGRF model gives a standard mathematical descrip- 
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Figure 3: The geomagnetic field lines, obtained using 
the IGRF(a), Tsyganenko 2001(b), and Tsyganenko 
2004 (c) geomagnetic field models. 
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Variation of plasma pressure with volume of the magnetic flux tube 




W, m/T 



Figure 4: Relationship between the plasma pressure, 
obtained from electron precipitating fluxes, measured 
by the Aureol-3 satellite and is the flux tube volume 
per unit magnetic flux obtained using the 1GRF (cy- 
cles), Tsyganenko 200f (squares), and Tsyganeneko 
2004 (diamonds) geomagnetic field models. 



tion of the Earth s main magnetic field and does not 
consider their distortion caused by magnetospheric 
current systems. Tsyganenko 2001 (TS01) and 2004 
(TS04) models represent the expected statistical re- 
sponse of the geomagnetic field to the orientation of 
the Earth s dipole axis, solar wind pressure, inter- 
planetary magnetic field, and appropriate geophysi- 
cal indices. TS04 model is a dynamical model of the 
storm-time geomagnetic field in the inner magneto- 
sphere that includes the temporal evolution of mag- 
netospheric current systems during the entire storm 
cycle. The input parameters are similar to those used 
by the TS01 model. Nevertheless a combination of 
the solar wind density, speed, and magnitude of the 
southward component of the IMF, averaged every five 
minutes, makes it possible to reproduce their evolu- 
tion during great geomagnetic storms. 

As it can be seen in Figure 3, the geomagnetic field 
lines differ significantly for all models, especially for 
the distances more than 4 Earth s Radii (Re). The 
TSOf model also has a singularity close to 6 Re in 



the tail direction. The same picture was observed 
for other two events, occurred at the end of main 
phase and in the beginning of recovery phases. For 
other events the singularities were observed at the 
distances more than 10 Re. Therefore, we analyzed 
the distribution of plasma pressure and other related 
parameters for the distances closer than the singular- 
ities often observed in TS01 and/or TS04 models. 

Stability of plasma pressure profile formed with re- 
spect to the development of interchange or flute in- 
stability can be an important factor that determines 
the processes during magnetic storms (see Antonova 
(2005, 2006) and references therein). It supposes, 
that increase of the plasma pressure inside the mag- 
netosphcrc in the region of low plasma parameter (3 
((3 — 2^op/£>2), where B is the value of the magnetic 
field in the equatorial plane, ?0 is the permeability 
of vacuum) by radial plasma transport and flux tube 
filling by ionospheric source is possible up to the limit 
determined by the interchange or flute instability de- 
velopment only. According to Kadomtsev (1963), the 
pressure profile should be lower that a critical limit 
given by . Here W — J dl/B is the flux tube volume 
per unit magnetic flux, dl is the element of field line 
length, B is the value of the magnetic field and in- 
tegration is done between conjugate hemispheres, 7 
is the ratio of specific heat and the value of const is 
determined by the boundary conditions. In case of 
adiabatic compression 7 = 5/3. Particle transport 
together with conservation of the number of parti- 
cles in the magnetic flux tube and first and second 
adiabatic invariants leads to nearly the same value 
7 = 7/4 (Tverskoy, 1997). 

The analysis of the stability of plasma pressure 
distribution in such a case requires the analysis of 
a relationship between the plasma pressure and the 
magnetic flux tube volume per unit flux P = P(W) 
(see Figure 4). The values of W were obtained by 
integrating the values of geomagnetic field along the 
field lines using the IGRF (circles), TS01 (squares) 
and TS04 (diamonds) models. As can be seen, there 
is a strong difference in the behavior of plasma pres- 
sure according to the IGRF and TS04 from one side 
and TSOf from other side. For the first two mod- 
els plasma pressure profile is situated at the dipole- 
like field lines, whereas TS01 has only three points 
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with negative pressure gradient at the dipolc-likc field 
lines. So we included 4 points from the tail- like con- 
figuration after singularity to be able to fit the pres- 
sure profile, as illustration. Nevertheless, we believe, 
that change from dipolc-like to tail-like configuration 
should be reflected in particle precipitations as well, 
what did not happen, indicating that the TS01 ge- 
omagnetic field configuration is much less probable. 
Radial pressure profiles were fitted by a power law 
P (W) W lw . Such fittings give the following val- 
ues for the exponents = —1-9 ± 0.1 for IGRF, 
lw = -0.21 ± 0.08 for TS01, and j w = -1.34 ± 0.07 
for TS04. It is easy to see, that the TS01 model is 
much more overstretched than the IGRF and TS04 
models. Nevertheless, it is reasonable to assume that 
at the distances less that 4Re, the difference between 
configuration of geomagnetic field lines and dipolar 
ones is not strong, and the correct value of 7 lies 
between -1.3 and -1.9. Unfortunately this accuracy 
is not enough to suggest specific plasma instabil- 
ity, which could be an important factor for the geo- 
magnetic storm development. Nevertheless, it allows 
tracking the relative changes in the plasma pressure 
profiles, related to the storm development. Figure 
5 shows the variation in the exponents for IGRF 
(black circles), TS01 (white squares), and TS04 (gray 
diamonds). As it can be seen, during February 26 
and 27 (days of year 57 and 58) IGRF model dis- 
plays strong decrease reaching unrealistic values of 
j w = —4.6. For this orbit, there are significant dif- 
ferences between IGRF from one side and TS01 and 
TS04 models from other side. This indicates that for 
this specific orbit (February 27, 1982 at 22 UT, ap- 
proximately) the geomagnetic field lines probably de- 
viated on dipolar configuration even at the distances 
of 6 Re. During the recovery phase all three mod- 
els indicate that the radial pressure profiles start to 
be flatter. All three models also indicate that at the 
end of main phase the maximum of plasma pressure 
was situated close to the Earth (L max = S.ORe for 
IGRF, Lmax=3.5 Re for TS01, L max = 3.2Re for 
TS04). It is interesting to mention that according to 
Tverskaya (1986), Tverskaya et al. (2003) lowest po- 
sition of the westward electro jet center Lmax during 
the storm is related to the maximum value of Dst- 



March 1-8, 1982 storm. Aureol-3 data. 
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Figure 5: From top to bottom: variation of Dst in- 
dex, variation of the jw, obtained due to fitting of 
P (W) W 7w , maximum pressure position in Re. Here 
black circles were used for IGRF, white squares for 
TS01, and gray diamonds for TS04 models 



variation Dst max in accordance with the relation: 



D st = 2.7510 4 L m 4 a:c nT 



(1) 



where Lmax corresponds to the position of the plasma 
pressure maximum. The relationship (1) was ex- 
plained by Tverskoy (1997), Antonova (2005, 2006) 
including the coefficient 2.7510 4 . For the storm ana- 
lyzed D st = 211 nT, that corresponds to the theoret- 
ical value L max — 3.34i?e. This value is in a rather 
good agreement with our estimation of the position 
of plasma pressure maximum despite we can analyze 
only part (may be only small part) of integral plasma 
pressure, excluding the contribution of more ener- 
getic particles. Interesting feature of the observed 
values of plasma pressure maxima appears when we 
compare these values with the solar wind dynamic 
pressure (Figure 6). It is possible to see a very good 
correlation, between the values of pressure maximum 
and solar wind dynamic pressure, excluding the main 
phase of magnetic bay, when the solar wind dynamic 
pressure is five times larger. We assume that this 
result must be tested in wider data sets, especially 
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Solar wind and plasma sheet pressure 
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Figure 6: Variation of the solar wind dynamic pres- 
sure from OMNI2 data base (black circles) and vari- 
ation of maximum plasma pressure in the inner mag- 
netosphere (white squares). 



considering great uncertainty of our methodology for 
obtaining values of plasma pressure maxima inside 
the magnetosphere. This result agrees with the re- 
sults of Borovsky at al. (1998) who established that 
the density of the plasma sheet is strongly correlated 
with the density of the solar wind. 

Correlation between the solar wind dynamic pres- 
sure and the maximum plasma pressure in the in- 
ner magnetosphere can be a consequence of magne- 
topause pressure balance during at least quiet peri- 
ods. Solar wind dynamic pressure is balanced by the 
pressure of the inner magnctosphcric magnetic field 
only in the regions of low plasma (3 parameter in- 
side the magnetosphere (for example, in the subsolar 
point). Nevertheless, the plasma [3 parameter can be 
higher in the near cusp regions and at the magneto- 
spheric flanks. Therefore, contribution of the inner 
magnetospheric plasma pressure in the total pressure 
balance in these regions can be very high. Taking into 
consideration, that in case of magnetostatic equilib- 
rium the plasma pressure is constant along current 
lines, it is possible to conclude that quiet time cur- 



rent lines that correspond to the region of plasma 
pressure maximum cross the magnetopause or close 
itself inside the magnetosphere in the near cusp re- 
gions, where plasma pressure must be close to the 
solar wind dynamic pressure. 

This result can be important for understanding 
the processes of the solar wind-magnetosphere cou- 
pling. Recently, Newell at al. (2006, 2007) pro- 
posed a "universal solar- wind magnetosphere cou- 
pling function" that represents the rate of change 
of the magnetic flux at the magnetopause and de- 
pends on the solar wind velocity, the value of in- 
terplanetary magnetic field, and the IMF clock an- 
gle: dF MF /dt = o 4 / 3 4 /3 «n 8/3 (a/2) . The func- 
tion proposed gives the best fits for all tested ground 
based and satellite indices, except the Dst. Neverthe- 
less, the authors argued "that using either the current 
value of p 1 / 2 or using p 1 ! 2 integrated over 72 hours 
improved Dst. ... this implies that only a portion of 
the advantage of "correcting" Dst for p 1 / 2 actually 
arises because the ring current perturbation on the 
ground-station magnetometers. In other words, the 
actual ring current itself does respond to p 1 ^ 2 , pre- 
sumably because of the increased proximity of the 
magnetopause" . Our results indicate about the pos- 
sibility of a direct influence of the solar wind dynamic 
pressure on the pressure in the inner magnetosphere. 

3 Discussion and Conclusions 

In this paper we present an analysis of the distri- 
bution of lower energy part of plasma pressure in- 
side the magnetosphere and reveal some characteris- 
tic features of plasma pressure evolution during ge- 
omagnetic storm. In particular, we found that the 
plasma pressure profile became steeper, especially 
near the end of the main phase, having the value near 
and even exceeding (in IGRF mapping) the maxi- 
mum steeping value 7n/=7/4 necessary for the de- 
velopment of interchange instability in the case of 
isotropic pressure and near to dipole magnetic con- 
figuration. Unfortunately empirical storm time mod- 
els of the magnetospheric magnetic field are not well 
developed. This fact is reflected in significant differ- 
ences in the plasma pressure profiles, obtained using 
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Tsyganenko 2001 and 2004 geomagnetic field mod- 
els, that demands additional studies of inner magne- 
tosphere pressure distribution using a self-consistent 
magnetospheric model in which the magnetic forces 
are equilibrated by plasma pressure forces, satisfying 
a condition of quasi-static equilibrium (such equilib- 
rium is realized when plasma velocity is much smaller 
than Alfven and sound velocities). One of the first 
important steps in such direction was made by Za- 
haria et al., (2005, 2006). Storm time plasma pres- 
sure profiles at Lj9 observed by AMPTE/SSE satel- 
lite (Lui et al., 1987) give the values of plasma pa- 
rameter (3 < 1. In our case all three used mod- 
els positioned the maximum of plasma pressure al- 
ways at quasi-dipole magnetic field lines, where the 
magnetic pressure was hundreds of times higher than 
the plasma pressure. We measure only a low energy 
portion of plasma pressure. Nevertheless, there still 
exists a probability of high plasma in the region of 
plasma pressure maximum. However, taking into ac- 
count the modeling results of Zaharia et al. (2005, 
2006) it is possible to suggest that low approxima- 
tion is valid at least for analyzed magnetic storm. 

It is interesting to mention that it was also found 
that the position of pressure maximum follows well 
the experimentally obtained dependence of the posi- 
tion of auroral electrojet on the value of Dst variation. 
This coincidence provides additional support of the 
theory of Dst formation during magnetic storms de- 
veloped by Tvcrskoy (1987), Antonova (2005, 2006). 

Interesting feature appears when we compare the 
values plasma pressure maxima with the solar wind 
dynamic pressure. A very good correlation is ob- 
served for all satellite passages, excluding the main 
phase of magnetic bay, when the solar wind dynamic 
pressure is five times larger. We plan to study this 
effect in detail in the nearest future. 
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